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Abstract 

We analyze the possibility of detecting one of the three scalar Higgs 
bosons in the next-to-minimal supersymmetric model at the future e + e~ 
linear collider, by examining their productions via the Higgsstrahlung pro- 
cess. The production cross sections of the three scalar Higgs bosons in 
e + e~ collisions are evaluated for the proposed cm. energies of the future 
e + e~ colliders, for the whole space of relevant parameters. We find that 
at least one of the three production cross sections is not smaller than 16 
fb, 4 fb, and 1 fb for yfe = 500 GeV, 1000 GeV, and 2000 GeV, respec- 
tively. Those numbers indicate that at least one of the three scalar Higgs 
bosons in the next-to-minimal supersymmetric model may be detected at 
the future e + e~ linear colliders via the Higgsstrahlung process. 



*E-mail: kim@physik.rwth-aachen.de (telephone +49-241-80-7239) 



The scalar Higgs bosons constitute an essential part of the standard model and 
the extended models therefrom [1]. To search for and discover them is thus one of 
the most important tasks of both theoretical and experimental particle physicists. 
As is well known, there is just one scalar Higgs boson, S, in the standard model, 
whereas there are multiples of them in the extended models. For example, the 
minimal supersymmetric standard model (MSSM) [2] has two scalar Higgs bosons, 
Si and S 2 , while the next-to-minimal supersymmetric standard model (NMSSM) 
[3] has three of them: Si, S 2 , and S 3 . 

The MSSM has just two guage doublet Higgs superfields, Hi and H 2 , and the 
NMSSM is a minimal extension of the MSSM, by introducing an additional gauge 
singlet Higgs superfield Af = (N, ip N , F N ) to the Higgs sector of the MSSM. Here, 
N is a Higgs singlet, ipN is a singlet higgsino, and Fn is an auxiliary field. As 
mentioned elsewhere [4], this singlet superfield can provide an economical way to 
avoid the so called //-parameter problem in the MSSM. Thus, while the MSSM 
introduces the dimensional /i-parameter by hand, the NMSSM can generate it in 
terms of the singlet superfield Af. 

In the NMSSM, the relevant superpotential is given by 

W = XHleH 2 N - \kAf 3 + ■■■ (1) 

where A and k are dimensionless parameters. We assume that they are real. The 
relevant soft-breaking part of the Higgs sector is accordingly given by 

V; oft = -\A x H^eH 2 N - hA k N 3 + h.c. + • • • (2) 

where A\ and A k are soft-breaking parameters having mass dimension. As the 
vacuum expectation value (N)q = x is developed, Ax emerges to correspond to 
the //-parameter in the MSSM superpotential W = fiHjeH 2 + — • • •. Note that 
the Peccei-Quinn symmetry is recovered if k = 0. Thus, AA; ^ 0. The two Higgs 
doublets Hi and H 2 give masses to the up-quark sector and the down-quark 
sector, when they develop the vacuum expectation values, Vi and v 2 , respectively. 

A side effect of the introduction of the new superfield is that it enriches the 
particle spectra of the NMSSM as well as introduces more parameters than the 
MSSM. While the charged sector of the NMSSM is the same as that of the 
MSSM, the neutral sector of the NMSSM has an extra neutralino, an extra scalar 
Higgs and an extra pseudoscalar Higgs boson. The number of parameters in the 
NMSSM is also larger than in the MSSM. For example, while one needs just two 
parameters in order to determine the Higgs boson masses in the MSSM, which 
can be chosen as /t and tan/5 = Vi/v 2 , six parameters are introduced for the same 
purpose in the NMSSM: A, x, k, A x , A k , and tan (3. 

The present experiments at LEP set lower bounds on the mass of the lightest 
scalar Higgs boson. The standard scalar Higgs boson may not be lighter than 65 
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GeV, while S\ in the MSSM may not have a mass below 45 GeV as it has not 
been detected yet [5]. For the scalar Higgs bosons in the NMSSM, it is shown 
that there are parameter regions where the lightest scalar Higgs boson Si may 
have a small mass without contradicting to the present experimental data [6]. 

Our main purpose in this article is to show that at least one of the scalar 
Higgs bosons in the NMSSM can be detected at the future e + e~ Linear Colliders 
with = 500, 1000, or 2000 GeV. This is done by calculating the minima of the 
production cross sections of the three scalar Higgs bosons in e + e~ collisions. We 
find that there are parameter regions in the NMSSM where at least one of the 
minima of the production cross sections are larger than the proposed discovery 
limit of the future e + e~ colliders, which implies that the future e + e _ colliders are 
capable of detecting at least one of the three scalar Higgs bosons in the NMSSM. 
A preliminary result of our analysis is reported elsewhere [7]. 

A remarkable result of the MSSM is that there is an upper bound on the mass 
of the lighter scalar Higgs boson given by 

msi < m z cos 2/3 

at the tree level. This is due to the fact that all quartic terms have gauge coupling 
constants. In case of the NMSSM with one Higgs singlet, there are quartic terms 
with couplings other than the gauge couplings. Thus, there is no simple upper 
bound on the mass of scalar Higgs boson as in the MSSM. The tree-level upper 
bound on m Sl is given as [8] 

m 2 Sl < (m£ ax ) 2 = m| (cos 2 2(3 + 2X " co f ° w sin 2 2 {3) (3) 

It approaches to the MSSM relation m 2 Sl < m| for tan/3 — > oo or A — > 0. 

The above relation may be cast into a simple form, by noticing that the 
coefficient of A above, 2cos 2 0vf/<?2> * s ver Y close to the neutrino- Z coupling 
constant g vv z = #2/2 cos Thus, denoting A = #2/ V2cos0w — V^g U uZ, 
one has \o(rriz) = 0.52 for g U vz(jnz) = 0.3714. With this value, the upper 
bounds on the Si mass at tree level is: m Sl < m z for A < A = 0.52 and 
m Sl < j^mz = 1.92 \m z for A > A [8]. 

The dependence of ms 1 on A comes from the fact that there is a quartic term 
with the coupling constant A in the NMSSM. It turns out that the upper bound 
on A may be relevant for ms 1 - An effective way of determining the upper bound 
on A is the renormalization group (RG) equation analysis [9,10]. 

The 1-loop RG equations involving A read as 

^ = JL(^ + 2A 2 3 2 _3 1 
dt 8vr 2V 2 * 2 y2 2 yi; 
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where t — In /j,, [i being here the renormalization scale. Numerical integration of 
the above RG equations by demanding the existence of no Landau poles up to 
the GUT scale yield the upper bound of A at the electroweak scale. We plot A max 
as function of tan/5 for some values of k, for m t = 175 GeV in Fig. 1 and for 
m t = 190 GeV in Fig. 2. 

We find that 0.64 < A max < 0.74 for 175 < m t (GeV) < 190. The upper 
bound on A is roughly independent of tan/3 for tan/3 > 3. Also, the two figures 
show that A max decreases with increasing k. The upper bound on k is about 0.7. 
Moreover, the lower bound on tan (3 is set at the electroweak scale. We have 
tan/3 > 1.24 for m t = 175 GeV and tan/3 > 2.6 for m t = 190 GeV. 

Using 0.64 < A max < 0.74, we obtain the tree-level upper bound on m Sl as 



Our results may well be compared with previous analyses. For example, the 
upper limit of A has been estimated to be about 0.87 [9] for h t > 0.5, where h t 
is the top Yukawa coupling, which yields ms 1 < 151 GeV. By similar analysis of 
RG equations, the upper limit of A was obtained as function of h t , thus yielding 
an upper bound on ms 1 as 140 GeV [10]. 

However, the above tree-level upper bound does not contain radiative cor- 
rections to the mass matrices of Higgs bosons. As in the case of the MSSM, 
the contributions from the radiative corrections is found to change the tree-level 
upper bound considerably. Several groups have calculated the higher-order con- 
tributions to the mass matrices of Higgs bosons and determined the corrected 
upper bound on m Sl [8,11,12]. 

We employ the effective potential method of Coleman and Weinberg [13] to 
evaluate the higher-order contributions at the 1-loop level. The 1-loop effective 
potential can then be written as [11] 



113 < m 



(GeV) < 131 



(5) 



V 1 (Q) = V (Q)+AV 1 (Q) 



(6) 
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where V (Q) is the tree-level potential evaluated with couplings renormalized at 
the scale Q, and 

1 A/f 2 S 

^W) = ^T r ^(l n — --) (7) 

is the 1-loop contribution. Here, M 2 is the field-dependent generalized mass 
matrix containing mass terms for all the particles in the NMSSM. 

The mass matrix M s of the scalar Higgs bosons at 1-loop level is, to a good 
approximation [11], given by the second derivative of Vi(Q) with respect to the 
Higgs fields. Taking only top quark and stop quark contributions into account, 
and assuming the degeneracy of the left and the right stop quark mass, we find 
that the higher-order contributions change the relevant elements of the tree-level 
mass matrix as [8] 



M\ 



12 



M, 



22 



M n + D n 

>12 



Mf 2 + D l 



M* 2 + D 
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where 



D 22 




(9) 



1 A\A\ 
6 M~ 



with At = —A t + Xx cot (3. The 6-quark contributions are not included as 
they are negligibly smaller than the t-quark contributions. Those higher-order 
contributions render the upper bound on m Sl as 

( m max )2 ^ (m max,l-loop )2 = (^2 + ^ ^2 q + ^ sin ^ + jj^ ^2 ^ (1Q) 



Incorporating the upper bound of the parameter A from the RGE analysis, we 
determine the upper bound on ms 1 in the parameter regions of the NMSSM for 
250 < x (GeV) < 1000, 250 < A t (GeV) <= 1000, and 250 < m~ t (GeV) < 1000, 
where x, A t , and are respectively the vev of the Higgs singlet, the soft breaking 
coefficient of top sector, and stop quark mass, and for 2 < tan/3 < 20, for 
175 < m t (GeV) < 190. Our result is: 

120 < m^ ax '^ loop (GeV) < 156 (11) 

for the lightest scalar Higgs boson in the NMSSM at the 1-loop level [8]. 
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Also, in terms of ra™ 8 *, it can be shown that the upper bounds on m S2 and 
ms 3 are given as 

( m max)2 _ Rj m ^ 

T^RJ 

( m ™x)2 _ {R 2 + R 2 )m 2^ 

i-(m+ Ri) 

where R\ = Uu cos (5 + Ui 2 sin f3 and R 2 = U 2 i cos (5 + U22 sin (5, Uij being the 3x3 
orthogonal matrix which diagonalizes the mass matrix M s of the scalar Higgs 
bosons. Clearly, R 1 and R 2 are complicated functions of the relevant parameters 
of the NMSSM. Note moreover that they satisfy the unitarity condition < 

R\ + R\ < l. 

That the upper bound on m Sl in the NMSSM at 1-loop level is between 120 
and 156 GeV suggests that the accessible region of the parameter space at LEP 
1 with yHi = mz might be very small. Actually, we have shown that the existing 
LEP 1 data do not exclude the existence of Si with m Sl = GeV [6]. 

For the future colliders with yfs = 500, 1000, or 2000 GeV, the situation is 
quite different. In this case, the production of the lightest scalar Higgs boson via 
the Higgsstrahlung process, e + e~ — > S±Z where Z decays further into a pair of 
fermion and antifermion, with real Z and real Si, is always possible, because the 
collider energy y/s is much larger than E T = m z + m Sl : 

212 < E™ ax = m z + m^ ax (GeV) < 248 

Thus, Et plays a kind of threshold energy and is an important quantity of our 
model. To obtain informations about how far the model could be tested, then, 
it would be helpful to derive a lower bound on the production cross section of 
the Higgsstrahlung process, which can be expressed as a function of the collider 
energy only. 

At the proposed center of mass energies of 500, 1000, or 2000 GeV for the 
future e + e~ linear colliders, the question is therefore not whether it is kinemati- 
cally possible to produce Si but whether the production rate is large enough for 
Si to be detected. If the production rate of Si be small, one should then examine 
as a next step if that of S2 or S3 is large enough. In order to by systematic, we 
consider all of the productions of Si, S 2 , and S 3 via the Higgsstrahlung process. 

The cross sections for the productions of the three scalar Higgs bosons via the 
Higgsstrahlung process can be expressed as 

01 (Ri , R 2 , m Sl ) = ctsm (msi ) R\ 

<T 2 (Ri,R2,ms 2 ) = a S M(m S2 )Rl (12) 
a 3 (R 1 ,R 2 ,m S3 ) = <J 5M (m 53 )(l - R\ - R\) 



ml < (m^ x ) 2 = 
m% < « 3 aX f = 
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where (JsmM is the cross section in the standard model for the production of the 
Higgs boson of mass m via the Higgsstrahlung process. Note that the production 
cross sections are also functions of ^/s of the collider, though not explicitly shown. 
A useful observation is that o"j(m™ ax ) < cr^mgj which allows one to derive the 
parameter- independent lower bound on <jj, as we will see shortly. 

Now, we calculate for given Ri and R 2 the cross section <Ji(Ri, R 2 ,m Sl ) for 
Si production via the Higgsstrahlung process. We repeat the calculations for S 2 
and S3 productions to obtain a 2 (R\, R2, ^s 2 ) and 03 (-Ri, R2, m s 3 )- Because the 
production cross sections decrease as the produced Higgs boson mass increases, 
one has 

a 2 (R 1 ,R 2: m S2 ) > a 2 (R u R 2 ,m^ x ) 

a 3 (R 1 ,R 2 ,ms 3 ) > <t 3 {R u R 2 , m^) (13) 

The second step is to calculate, by allowing ms 1 to vary from its minimum 
to its maximum value, for given Ri and R 2 , the set of the three production 
cross sections <7i(i?i, R 2 , m Sl ), a 2 (Ri, R 2) m™ ax ), and a 3 (R 1: R 2: m™ ax ). Recall 
that both m™ x and m™ ax are functions of m Sl . 

It is quite clear that each of the three production cross sections will exhibit its 
own minimum (and maximum) for a certain value of ms 1 in the range between 
its lower and upper bound. Naturally, ai(R\, R 2 , mgj is minimized when ms 1 is 
maximized. That is, the minimum of a\(R\, R 2 , rngj is given by a\(R\, R 2 , m™ 8 *). 
Let us denote it by ^(.Ri, R 2 ). 

For a 2 (Ri, R 2 , m™ 2 ax ) and a 3 (Ri, R 2 , m™^), it is not straightforward to tell 
where their minima occur within the allowed range of mg 1 , because both m™ ax 
and m™ ax are dependent not only on m Sl but also R\ and R 2 . Nontheless, the 
minima of o 2 (R\, R 2) m™ x ) and a 3 (i?i, R 2: m^) would exist for certain values of 
m Sl . Let them be denoted by cr™ m (i?i, R 2 ) and a 3 nm (Ri, R 2 ), respectively. (Note 
that the value of ms 1 that yields <7™ in (i?i, R 2 ) is generally different from the one 
that does a 3 Tim (R 1: R 2 ), for given Ri and R 2 .) 

Then, by comparing the three minima with each other, we can establish 

a(R u R 2 ) = max{ < 7 1 min (i2i, R2), < n (Ri, Ri), < in (^i, ^2)} (14) 

for given Ri and R 2 . In other words, 

<n(Ri,R2,m Si ) > <T i (R 1 ,R 2 ,mH* K ) > af*(R u R 2 ,m St ) = a(R,,R 2 ) 

for some %. The meaning of &(R±, R 2 ) cannot be missed: At least one of the three 
scalar Higgs bosons has its production cross section via the Higgsstrahlung process 
in e + e~ collisions larger than a(Ri,R 2 ) for a certain set of relevant parameters 
of the NMSSM. 
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The third step is to plot a(R±,R 2 ) in the i? 2 )-plane for given y/s. The 
triangular area in the plane defined by < R\ < 1, < R\ < 1, and < 
(1 — R\ — R\) < 1 represents the whole space of the relevant parameters of the 
NMSSM, because it is the entire physical area in the (Ri, -R2)-plane. 

We find that a(Ri,R 2 ) vanishes for the collider energy y/s less than E T = 
m max _|_ m ^ This is the case for LEP II with yfs < 205 GeV, which implies that 
even the highest proposed energy of LEP II would most probably not be able to 
test our model completely. If \fs is larger than Et, then a(Ri, R 2 ) never vanishes 
for that collider energy. Thus, for given y/s > E T , a(Ri,R 2 ) exhibits a nonzero 
minimum somewhere in the (Ri, R 2 )-pl&ne. 

The final step is to establish the minimum of cr{R\, R 2 ) in the i?2)-plane 
for given \fs. Let it be denoted as <?o(y/s). Thus, we have 

a{R u R 2 ) > a (^7s) 

This minimum is a universal minimum in the sense that it is independent of 
the parameters of the NMSSM. It is the absolute lower bound on at least one of 
(Ji(Ri, R 2 , msj > o"o(v^) f° r some i: At least one of the three scalar Higgs bosons 
may be produced via the Higgsstrahlung process in e + e~ collisions with produc- 
tion cross section larger than o"o(a/s). This minimum is therefore a characteristic 
quantity of the NMSSM. 

Let us now plot a(R\,R 2 ) in the (R\, i?2)-plane for given \fs in order to 
evaluate a^y^s), its minimum. For simplicity, we set mc? ax = 145 GeV. The 
dependence of a(R 1 , R 2 ), and hence a (y/s), on m<? ax is quite small for y/s > 500 
GeV, as will be shown shortly. 

In Fig. 3, we plot cr(Ri, R 2 ) for y/s = 500 GeV. The minimum is found to 
be about 16 fb. When the discovery limit is about 30 events, one would need a 
luminosity of about 25 fb for the future e + e" collider, which is a realistic one. In 
Fig. 4 and Fig. 5, we repeat our plottings for ^/s = 1000 GeV and 2000 GeV, 
respectively. We find that the minimum of the production cross section is 4 fb 
for v/i = 1000 GeV and 1 fb for 2000 GeV. 

To see the dependence of ^(y^) on the collider energy, we plot it as a function 
of t/s in Fig. 6 for some values of m™ 8 *. It is observed that <Jo(\fs) becomes 
largest for around y/s = 300 GeV. Moreover, the effect of m™ x is most dominant 
there. For y/s > 500 GeV, a (y/s) decreases as y/s increases. Also, the effect of 
m™ on cr (y/s) becomes rather small there. 

As an illustration, the individual production cross sections <Jj(i?i, R 2 , m Si ) 
are calculated for an exemplary set of the relevant parameters of the NMSSM. 
Here, not only the Higgsstrahlung process e + e~ — > S%Z — > Sibb but also the 
contributions from two other important processes are considered. The two other 
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processes are the process e + e~ — > Z — > 66 — > S^frfr where 5, is radiated off from 
6 or 6, and the process e + e~ — > Z — > S^P, — > S'jfefe where P, (j = 1,2) is a 
pseudoscafar Higgs boson. 

The resufts are plotted in Fig. 7 and Fig. 8, as functions of the collider energy 
y/s, for A x = 220 GeV, A k = 160 GeV, x = 1000 GeV, A = 0.12, k = 0.04, and 
tan j3 — 2. The production cross sections in Fig. 7 are the tree-level ones, whereas 
those in Fig. 8 are those with the 1-loop corrections via the effective potential. 
The relevant parameters for the 1-loop corrections are taken as A t = 0, = 1000 
GeV, and m t = 175 = GeV. 

For this exemplary set of parameters, the 1-loop corrections is seen to be 
important for the collider energy of about 150 GeV and decreases as the collider 
energy increases. Also, for y/s = 500 GeV and beyond, we find that S 2 would 
most dominantly be produced for the exemplary set of parameters, both at the 
tree level and the 1-loop level. 

The possibility of detecting one of the three scalar Higgs bosons in the NMSSM 
at the future e + e~ linear collider has been analyzed by examining their produc- 
tions via the Higgsstrahlung process. The whole space of the relevant parameters 
of the NMSSM is conveniently represented by a triangular area in the (Pi,Pl2)- 
plane, where the production cross sections of the three scalar Higgs bosons in 
e + e~ collisions are evaluated for given cm. energy of the e + e" collider. 

For the three production cross sections, we have searched their minima in the 
i?2)-plane. Then, from the three minima of the production cross sections, 
we have established a characteristic quantity for them, a (y/s), which is a kind 
of universal minimum of the production cross section. Actually, it is the largest 
one among the three minima of the production cross sections. We have obtained 
that a (y/s) = 16 fb, 4 fb, and 1 fb for = 500 GeV, 1000 GeV, and 2000 GeV, 
respectively. As stressed in the preceding section, these numbers indicate that at 
least one of the three scalar Higgs bosons may be produced via the Higgsstrahlung 
process in e + e _ collisions with production cross section larger than cr (y/s). 

What we have found indicate that, for a certain set or regions of the relevant 
parameters of the NMSSM, at least one of the three scalar Higgs bosons might 
have a reasonable mass and a reasonable production cross sections such that it 
might be produced at the future e + e~ colliders. We conclude that the Higgs sector 
of the NMSSM can most probably be tested at the future linear e + e~ colliders 
with v/i = 500, 1000, or 2000 GeV. 

After we completed this article, we have been informed that a similar conclu- 
sion has been drawn by Gunion, Haber, and Moroi [14]. However, their approach 
is quite different from ours. 
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Figure Captions 



Figure 1: The upper bound on A as function of tan/3 for some values of k, for 
m t = 175 GeV. 

Figure 2: The upper bound on A as function of tan/3 for some values of k, for 
m t = 190 GeV. 

Figure 3: The contours of a(R 1 ,R 2 ) in the (R u i? 2 )-plane for m^ x = 145 GeV 
and v/i = 500 GeV. 

Figure 4: The contours of a(Ri, R2) in the _R 2 )-plane for m™^* = 145 GeV 
and v/i = 1000 GeV. 

Figure 5: The contours of a(R 1 ,R 2 ) in the (R 1: i? 2 )-plane for m™ x = 145 GeV 
and v/s = 2000 GeV. 

Figure 6: <J (y/s), the minimum of <r(Ri, i? 2 ), as function of the collider energy 
y/s for some values of m™ x . 

Figure 7: The production cross sections at the tree level <r(e + e~ — > bbSi) (i = 
1, 2, 3) as function of the collider energy ^/s, for A\ = 220 GeV, Ak = 160 GeV, 
x = 1000 GeV, A = 0.12, k = 0.04, and tan/3 = 2. The solid, dashed, and 
dashed-dotted curves correspond respectively to the Si, S2, and 5*3 productions. 

Figure 8: The production cross sections at the 1-loop level <r(e + e~ — > bbSi) 
(i = 1, 2, 3) as function of the collider energy a/s, for A\ = 220 GeV, A^ = 160 
GeV, x = 1000 GeV, A = 0.12, k = 0.04, and tan/3 = 2. The 1-loop parameters 
are taken as A t = 0, m t ~ = 1000 GeV, and m t = 175 GeV. The solid, dashed, and 
dashed-dotted curves correspond respectively to the Si, S 2 , and S3 productions. 
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Figure 6: ao(y/s), the minimum of cr(i?i,i? 2 ), as function of the collider energy 
y/s for some values of m™ ax . 
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Figure 7: The production cross sections at the tree level <r(e + e~ — > bbSi) (i = 
1,2,3) as function of the collider energy y/s, for A\ = 220 GeV, Ak = 160 GeV, 
x = 1000 GeV, A = 0.12, k = 0.04, and tan/3 = 2. The solid, dashed, and 
dashed-dotted curves correspond respectively to the S\ } S2, and S3 productions. 
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Figure 8: The production cross sections at the 1-loop level <r(e + e~ — > bbSi) 
(i = 1,2,3) as function of the collider energy y/s, for A\ = 220 GeV, Ak = 160 
GeV, x = 1000 GeV, A = 0.12, k = 0.04, and tan/3 = 2. The 1-loop parameters 
are taken as A t = 0, = 1000 GeV, and m t = 175 GeV. The solid, dashed, and 
dashed-dotted curves correspond respectively to the S\ } S2, and S3 productions. 
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